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Abstract 
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1 Introduction 

The goal of this paper is to extend to the g-discrete orthogonal polynomials 
on the non-uniform lattice x(s) = (q s — 1) / (q — 1), the characterization for 
classical orthogonal polynomials in terms of the standard Stieltjes function (see 
|23j for original and historical issues on the Stieltjes function; for modern and 
advanced notions see [8], [15], [25] and the references therein). 

The properties of q-Charlier, q-Kravchuk, q-Meixner and q-Hahn orthogonal 
families are well known (see [H[5D] and [H]): they satisfy a hypergeometric- type 
difference equation pTJ]) ; their finite differences constitute an orthogonal poly- 
nomial family ([T5|; they can be expressed by a Rodrigues-type formula (fT5|) ; 
their associated orthogonalizing weights satisfy a Pearson-type difference equa- 
tion (Tl~3]l - (ri4l) ; they verify a three-term recurrence relation (fT7|) . These proper- 
ties, among others, characterize the classical orthogonal polynomials (see e.g. 
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[TJ[T21[T11[T7])- Most of these characterization properties have been already ex- 
tended to the g-polynomials. However, special attention must be given to [19], 
particularly for the unified study of some classical ^-polynomials carried out. 
An algebraic approach developed by Maroni [T^ is used as a crucial ingredient. 
Indeed, the starting point is a distributional difference equation fulfilled by the 
moment functional with respect to which the aforementioned ^-polynomials are 
orthogonal. The motivation for studying the problem considered in the present 
paper comes precisely from an unproven assertion given in Proposition 2.27 (b) 
of the aforementioned paper [15] . In this paper we will prove the above men- 
tioned assertion for classical g-polynomials on the lattice (q s — l)/(q— 1), i.e., a 
(/-moment functional IA is classic if and only if the Sticltjcs function associated 
with it, given in terms of the q- falling factorial, solves a non- homogeneous ver- 
sion of the Pearson-type difference equation associated with U (see Theorem [T] 
below). 

The advantages of the proposed approach over existing works consist of 
working directly with linear functionals rather than power series expansion for 
the Stieltjes function, i.e., of looking for the intrinsic properties that may satisfy 
the linear functional under consideration acting, either on the vector space of 
formal power series or on the vector space of polynomials. Our characterization 
provides a general framework which can be extended -like in the continuous 
case- to the study of g-semiclassical orthogonal polynomials [T5]. In addition, 
for illustrating our main result, the Stieltjes function given in terms of the q- 
falling factorial for the four canonical cases, namely, q-C'harlier, q-Kravchuk, 
q-Meixner, and q-Hahn is explicitly calculated and expressed in terms of hyper- 
geometric functions. As far as we know, there were not explicit hypergeometric 
expressions for these Stieltjes functions in the literature. 

The structure of the paper is as follows. In Section [5J we give a quick 
overview of the relationship between orthogonal polynomials and Stieltjes func- 
tion. In Section [3] we compress some necessary definitions and tools. Lastly, 
in Section 2] the Stieltjes function for g-orthogonal polynomials in terms of the 
q -falling factorial basis is introduced and the main theorem for the proposed 
characterization is proved. This result constitutes a (/-analogue of the characte- 
rization given in [TS] for classical discrete orthogonal polynomials (see also [S]). 
Finally, in Section [5] the explicit expressions for the corresponding solutions of 
the difference equation that characterizes the classical g-orthogonal polynomials 
on the non-uniform lattice x(s) = (q s — l)/(q — I), i.e. the Stieltjes functions, 
are studied. 

2 Stieltjes function and orthogonal polynomials 

We start by reviewing some results needed for the sequel; mainly extracted from 
EOl [22] . The linear space of polynomials with coefficients in C (the set of 
complex numbers) is denoted by P, and its dual space by P*, whose elements 
are linear functionals (moment functionals). 

Let the action of U € P* on p € P be denoted by the duality bracket (U,p). 
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The moments of U are given by (U, x k ) — Uk, k = 0, 1, . . ., and by linearity yields 
the following relation (U,P n (x)) = J2k=o a k u k, where P n (z) = YX=o a k zk - 

Formally, if U is a moment functional determined by the moment sequence 
{uk}k>o, the Stieltjes function associated with U can be defined as the gener- 
ating function of these moments by means of the series expansion 

fc>0 

Notice that S(z) corresponds to its sum when the Laurent expansion (TIJ con- 
verges for some z; otherwise S(z) represents its analytic continuation -provided 
that it exists-. 

The above function ([1]) can be approximated by rational functions with pre- 
scribed order near infinity |22) . This approximation problem -often known as 
(n — 1, n) Pade approximation near infinity- consists in finding the polynomials 

Pn{z) = Ylk=o a k zk and Qn-i(z) = J2k=o Pkz k , such that 

P n (z)S{z)-Q n - 1 (z) = 0{z- n - 1 ) 1 (2) 

where degP„ < n and degQ„-i < n — 1. 

By getting Q n -i{z) equal to the polynomial part of P n (z)S(z) and using 
the prescribed order at infinity in the interpolation condition ([2]). the following 
linear system of n homogeneous equations [22] 

n 

Uj+jCii = 0, j = 0, . . . ,n - 1, (3) 

i=Q 

for determining the n + 1 unknown coefficients of P n (z), holds. Hence, P n (z) is 
determined up to a multiplicative factor. Observe that the linear system is 
equivalent to the following orthogonality conditions: 

(U 1 P n {x)x k ) = 0, fc = 0,...,ra-l. (4) 

Consequently, P n (z) is orthogonal with respect to any polynomial of degree less 
than n for the linear functional U. On the other hand, we have selected Q n -i(z) 
equals to the polynomial part of P n (z)S(z), hence one has 

o.-.m=(m /- ( ']:* w ). (5> 

where it is assumed that U acts on the variable x. Notice that ([S]) leads to the 
formal relation 

P„0) (U, 7— - Qn-l(z) = (U, ( 6 ) 

\ (z — x) / \ z — x I 

which indeed implies -at least formal- the linearity extension of the linear func- 
tional U to the vector space of formal power series. Moreover, regardless of the 
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convergence of ([I} , it is straightforward to see that this series provides a means 
of approximating (hi, (z — Indeed, for every nonnegative integer n, 

(".^i)-(".|^) + <«.(j)"^i)- 

Thus, for z — > 00, in every sector e < argz < n — e (0 < e < ir/2) one gets 
\z — x\ > \z\ sine. Hence, 

(«.^i)-E^r-0.(^). 

N ' fe=0 

which holds for all z in the above sector. Analogously, for the remainder approx- 
imation term we get the prescribed order z~ n ~ 1 near oo, in the above sector, 
i.e 

\ z — x I 

Moreover, we recall that in the framework of the algebraic approach devel- 
oped by Maroni [THl H3 [TB] it is established a topological isomorphism between 
the space P* and the space of formal power series endowed with appropriate 
topologies. Indeed, it is always possible to associate (P n )n>o of degree n with 
a unique sequence (U n ) n >Ot U n £ P*, called the dual sequence of (P n )n>o, such 
that (U n ,P m ) = 0, for n ^ m, and equals 1 for n = m. Indeed, any functional 
WeP* can be represented as follows 

n>0 

which makes use of the linearly extension of U to the vector space of formal 
power series. Therefore, the Stieltjes function associated with U in the Pade 
approximation problem j2j has the formal representation 

S{z) = (U^z-x)- 1 ) . (7) 

From now on we will consider a specific class of moment functionals. Let fj, 
be a Borel measure on the real line R (with infinitely many points of increase) , 
supported on Q C M. We define the linear functional on which we will focus our 
attention 

(U, P ) = [ P (x)dfi(x), P e¥. (8) 

In particular, we will consider 

N 

H = ^2p(x)S Xk , p(x k ) > 0, Xk G K, iVeNU{+oo}. 

k=0 

This measure p is discrete (with finite moments u k — J n x k dfi(x), k > 0) 
formed by a linear combination of Dirac measures at the points xq, . . . ,xn- By 
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N we denotes the set of all nonnegative integers. Therefore, the orthogonality 
condition ([4]) for the ^-polynomials on the lattice {x(s) i-> R + : s = 0, . . . , N} 
with respect to §8§ is defined as follows -see [2H e 1- (4-7), p. 256] as well as [2"T] 



N 

(U,P n x k )=Y / Pn(x(s))x k ( S )p( S )Ax(s-l/2) = 0, fc = 0,...,n-l, (9) 

where, Ai(s) = x(s + 1) — x(s) denotes the forward difference operator. In 
general, the polynomial P n (x(s)) is called discrete orthogonal polynomial. In 
what follows we will denote any polynomial P n (x(s)) simply as P„(s). 

The term 'classical discrete orthogonal polynomial' [20j (also known as clas- 
sical orthogonal polynomials of a discrete variable) obeys the fact that x(s) = 
c\q s + C2q~ s + c 3 , (q € M + \{1}) or x(s) — c 4 s 2 + c 5 s + c 6 , where the constants 
Ci G M (i = 1, . . . , 6) are independents of variable s, and the orthogonalizing 
weights p(s) are solutions of the Pearson-type difference equation (see formulas 
JT3l)-(IT4j) below and [2QI pp. 70-72]). 

A remarkable feature of the classical discrete orthogonal polynomials is that 
they satisfy the hypergeometric-type difference equation 

a{s) Ax( s -l/2) V x(s) + T{3) AxV) + Ay(s) = °» (10) 
a(s) = a(x(s)) - (f(x(s)) A x(s - l/2))/2, t(s) = f(x(s)), 

where \/y(s) = Ay(s — 1) is the backward difference operator, and 

a(s) = a 2 x 2 (s) + axx (s) + a (11) 
T{s) = b lX (s) +b Q . (12) 

Equation (jlpp is a discretization of the very well known hypergeometric differ- 
ential equation (see [20] ) 

a(x)y"(x) +t(x) + Xy(x) = 0, with dega < 2, degf = l, kl. 

Notice that the solution p(s) of the Person-type equation (written in two equiv- 
alent forms) 

A[a(s)p(s)]=r(s)p(s), A= . . A 1/9V (13) 

Ax(s — 1/2) 

V[(a(s) + r(s) V x(s + l/2))p(s)] = r(s)p(s), V d = f Vx J +1/2 y (14) 



allows to write (flO)) in self-adjoint form. This symmetrization factor of (| 10[) 
with the condition a{s)p{s)x k (s — l/2)| s _ Q N — 0, fe = 0, 1, . . ., is called the 
orthogonalizing weight of ©• This boundary condition is compatible with the 
existence of all moments of ([9]). 

It is well-known that many properties of orthogonal polynomials can be 
deduced from (|9|); in particular, equation (|10[) . In turn from this equation 
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follow important properties like the Rodrigues-type formula (see [30] [H] for 
general theory and [2] for main data of g-polynomials on non-uniform lattice 
x(s) = (q s -l)/(g-l)). 

P n (s) = ^-vW[p n (s)}, VW = V , , V , , (15) 

p(s) Lr V^l(s) V^(s) VZn(s) 

where a;fe(s) = x(s+ |), and /o„(s) = p(s + n) Yim=i c(s + m); the orthogonality 
of the differences 



N-l 



> — : — a — t-^Pi(s) A x(s) = n 

t-~t At s Ar ^ W W 



s=0 



As(s) Aat(s) 



AP n ( s ) 



Ax(s) 



2 



(16) 



where ||-|| denotes the norm in the Hilbert space L^; the three-term recurrence 
relation 

x(s)P n (s) = a„P„+i(s) + /3„P n (s) + 7„P n _ 1 (s), n = 0, 1, . . . ,N, (17) 

where P_i — 0, Po is a real constant, and the coefficients a n , /3„, and 7„ are 
given explicitly in terms of (TTT1) and (TT2"]) (see [2], expressions in (42)). 

A second linearly independent solution of equation (1101) -for a specific choice 
A = A„, see [20]- is the so-called function of the second kind on non- uniform 
lattice (see [24], formula (5.1)). This function is also a second linearly indepen- 
dent solution of equation (fT7|) . For a relationship between the function of the 
second kind and the Stieltjes function we refer to [20] [24] . 



3 Basic definitions and notations 

Here, we will use the g-analogue of the Pochhammer symbol [TT] [T2] [5D] 

{a;q) k = JJ (l-aq j ), for k > 0, and (a;q) = l. (18) 
o<j<fc-i 

By Pochhammer symbol (s)k and falling factorial [s]k we mean 

(*) fc = s(«+l)---(* + fc-l), (s)o = l, k>l, 
[s] k = (-!)*(-*)*. 

respectively. Notice that the subscript in the above expressions is a nonnega- 
tive integer, whereas in the sequel symbol [s] q denotes the g-number defined as 
follows 

For the g-factorial we use the following definition 



! — 



< [n -!}„■■■ [2]„ [1]. , neN, [01! = 1. 
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Observe that in [111 p. 7] the g-number is defined in a different way. 
The q-hypergeometric series r ip s is defined as 



/ ai,. 


. . , a r 










V h,- 


..,b s 





\ - («i;g) fc ---K;g) fc r n * (*) 



{h;q) k ---{b s ;q) k L - (g;gV 



1 + s-r 



where (™j denotes the binomial coefficient. 

In what follows, the g-falling factorial defined as 



JJ x{s-j), for fc>0, and [s} { °> = 1, (19) 

0<j<fc-l 



where x(s) — , will play the same role that the power x k plays in a formal 

series expansion. Indeed, the q-falling factorial basis {1, [s]^ 1 \. . . , [s]^ , ... } 
has been suggested to be more natural than the basis {1, x(s), . . . , x n (s), . . . }, 
mainly when multiple orthogonal polynomials are considered [5]. Observe that 
the q -falling factorial [s]^ is a polynomial of degree k in x(s). Accordingly, 
the action on it of the forward and backward difference operators defined in 
d-d gives 

A [.]« = q^~ k [*]« , and V [„]« = ^ [*]« [ a - 1]^ , 

respectively. 

In addition, the g-falling factorial can be rewritten in terms of the q -analogue 
of the Pochhammer symbol as follows 



(fe)= (g s ;g) fcgfc ( s -*=i) , , ri(0) = 



(<z-i) fe 



q k \ s - — ), k>l, and [s]^ = l. (20) 



Clearly, the Pochhammer symbol and the falling factorial can be recovered as a 
limiting case 

lim to?' = (-l)*H0* = [*]*■ 

q— >1 H 

Let § be the spanning set of {1, [s]^ , . . . , [s]^ , . . . } over R. Observe that § 
coincides with the linear space of polynomials of discrete variable x(s) with real 
coefficients. §„ denotes the corresponding linear subspace of dimension n. Let 
us define the algebraic dual space of S, denoted by S* , as 

S* = {U : § — > R, such that U is a linear functional}. 

The real number it* = (lA, [s\ q n ^, n > 0, is said to be the ^-moment of order 

n of U and the sequence {u* } n >o is called the g-moment sequence associated 
with U. 
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Definition 1 Let {u-?}„>o be a sequence in M. IA is said to be the q-moment 

functional associated with the sequence {u^} n >o, if (IA, [s]^) = u^ L and it can 
be extended to § by linearity, i.e., if 

n n 

r(s) = £ofc[s]W, then (Z/,r) = ^ o*«g. 
Accordingly, from Definition [T] and equation one has 

N 

(M,r)^ r ( S )p( S )Ai( S -i), (21) 
where, the g-moment of order k is given by 

N 

< = (U, [s}f ) = £ Pis) A x( S - i), fc = 0, 1, . . . (22) 

s=0 

Below we sketch some definitions and results that we will use in the next sections 
(for comprehensive definitions and theorems on g-moment functionals we refer 
to [19]). 

Any moment functional U, which maps S into § has a transpose V T : S* i— > 
S* defined by (V T U,p) = (U, Vp), Mp G §, MIA G §*; then AU, VIA, and are 
defined by duality according to the following definitions. 

Definition 2 Let U be a q-moment functional. The forward and backward dif- 
ferences of the moment functional IA, denoted by AlA, and VIA, respectively, are 
defined as follows 

(AU,p(s)) = -(U,Vp(s)), where p e §, 
(VU,p(s)) = -(lA,Ap(s)). 

Here, the forward difference operator A on moment functionals is minus the 
transpose of the backward difference operator on polynomials. 

Definition 3 Let p G §. The linear functional V = plA is said to be the 
left-multiplication oflA by a polynomial p if 

(V,r) = (U,pr) , reS. 

Definition 4 // (cr, r) are polynomials of minimum degree such that 

A(o-IA) = tU, (23) 

define the class of IA as the nonnegative integer number s such that 

s = maxjdegcr — 2,degr — 1}. 

The polynomial sequence {Pn(s)} n>0 orthogonal with respect to the q-moment 
functional IA of class zero is said to be classical orthogonal polynomial sequence. 

Equation ([23]) is a q- analogue of the functional Person equation [TUJ [15] . 
In general, the polynomial sequences orthogonal with respect to such moment 
functionals IA are said to be semiclassical orthogonal polynomial sequences. 
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4 Stieltjes function for g-orthogonal polynomi- 
als in terms of the g-falling factorials 

Assume that IA is the moment functional determined by the q -moment sequence 
{ u k}k>o- Below, in the same fashion that the continuous case ([T]) we will define 
a q-analogue of the Stieltjes function associated with U (in short g-Stieltjes 
function) as the generating function of these g-moments by means of a series 
expansion, but using the g-falling factorial basis (jT9j) instead. More precisely, 



Definition 5 LetlA be a q-moment functional defined on the linear space S. The 
q-Stieltjes function associated with U is defined as the formal series expansion 



k>Q 1 



The following result deals with the (/-analogue of the formal representation 
([7]) . Indeed, we will use a q-analogue of the Chu-Vandermonde convolution (see 
[12] for details) 



2<Pl 



q; q n c/b 



(cb x ;q) n 



(c; q) n 

In particular, one can select s € {0,1, ... , N} and z £ C such that 



2fl 



q s ,q 



q l-z 



We will use this expression in the proof of the following lemma as well as the 
linearity extension of ti to the linear space of formal power series in accordance 
with the discussion given in section [2] with the usual adaptation required for 
discrete measures. Indeed, using the relation (see p. 262 [2~i] ) 

1 _ J_ V ( x ^\ k . x "( s ) 
77) 2^ 



x(s) — x(z) x(z) ^— J \x{z) J [x{s) — x(z)]x n (z) ' 



and proceeding like in ([6]) one gets an asymptotic expansion of (|30[) as x(z) — > oo. 
However, we refer to [161 1171 118) for a more comprehensive analysis. 



Lemma 1 The q-Stieltjes function \2J$ admits the following formal represen- 
tation 

S q (z) = (u,[x(z)-x( S )}- 1 ), (26) 
where it is assumed that ti acts on the discrete variable x(s). 
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Proof. Relations (f2T)|) and (j2"2")l enable us to rewrite (fM|) as follows 

_ (u, [*]< fc) ) _ (q - 1) («, (g- s ;g) fc ^-(")) 

Particularizing the well known relation (a; q) k+n = (a; q) n (aq n ; q) k -see [TT]- for 
the choice a = q~ z and n = 1, i.e. 

(<T z ;<z) fe+1 = (g 1- *;?)^ 

it follows that equation (|27p can be rewritten as 



„ , g-1 ^ {U,q k ^Hq~ s ;q) k ) 



If one uses (|2"5f equation (|2"5)) transforms into the expression 

*«-("-i5)&3&>- (29) 

Taking into account relation (jT8j) , a straightforward calculation leads to the 
following equalities 

1 (q- z ;q) s (l-q)(l-q-*) 1 



x(z) 



x(z){q 1 - z ;q) s (1 - <f)(l - q-*+ s ) x{z) - x(s) ' 
which makes expression (|2^|) coincident with (|26p ; therefore the statement holds. 
Remark 1 From 1121]) and 126\) one can write \2m up to a constant factor as 

S , W .f?™, (30 , 

s=0 w w 

which reveals its singularities. 



The Stieltjes function (|30j) is related to the function of the second kind 
on non-uniform lattice (see [24], formula (5.1)) for the specific choice n = 0. 
This function of the second kind is the second linearly independent solution of 
equations (fT0|) and (fT7|) . respectively. 

In the next results the subscripts are used to indicate on which variable the 
difference operators are acting. 

Lemma 2 The following relation 

s ' x(z) — x(s) / \ ' z x(z) — x(s) / ' ^ ^ 

holds. 
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Proof. By using the operator V defined in (fT3)| one can compute the back- 
ward difference of the Stieltjes function given in the q-falling factorial basis. 
Indeed, by Lemma Q] one gets 



VSq{z) = V<z + h) [Sqiz) ~ ^ ~ 1)1 



1 



or equivalently, 



x(z) — x(s) I \ x{z — 1) — x(s) 



VS q (z) = {U,V z ' 



' x(z) — x(s) 

A straightforward computation leads to the following relation 

1 A 1 

V z = — = A s 

x(s) — x(z) tt(s,z) x{z) — x{s) ' 

where 

7r(s, z) = [x (s + 1) — x (z)] [x (s) — x (z)] 

= q[x(s)-x(z-l)][x(s)-x(z)], (32) 

is a polynomial of degree two both in the variable s and z. 
Hence, 

vs ^ = -( u ^^zTWj)' 

which implies that equation (|3"Tj) is fulfilled. 

The following theorem establishes the main result in this paper. 

Theorem 1 The polynomial sequence {P n {s)}n>o orthogonal with respect to 
the moment functional U is classical if and only if the Stieltjes function \2J$ ( in 
terms of q- falling factorials) satisfies the following first order non-homogeneous 
ference equation 



V[(a(s)+T(s)vx(s + l/2))S q (s)]=T(s)S q (s) + C q , C q eR\{0}, (33) 
where the constant 

C q = (a 2 <T 1/2 + \biq-\q - 1) - fci) ug, 

depends on the polynomial coefficients a 2 and b\ ofa(s) and r(s), respectively 
(see formulas JTTp -IXI j) ). 

Proof. Let us introduce the linear functional V = p{s)lA, where 

P ( s ) = a ( s ) +t(s)\7 x(s + i) 
= ax 2 (s) + fix (s) + 7, 
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and 



a = a 2 + ~hq 1/2 (q - 1) , 

P = at + \biq- 1 ' 2 + \b q- 1/2 (9-1), 

7 = a + ^b a q~ 1/2 . 



Here we have used that 

\jx{s + ^)=q- 1 ' 2 (g-1) 



z(s) 



1 



9-1 

=g-V 2 {q-l)x{s)+q-V 2 . 



Notice that 



x (s) — x (s — 1) 
a;(s+l/2)-a;(s-l/2) 



-1/2 



Hence, 

Vp (s) = aVx 2 (s) + f3Vx (s) 



a 



x 2 (s) - x 2 (s - 1) 



x (s + 1/2) - x (s - 1/2) + 
a<T 1/2 [as (s) + x (s - 1)] + ^" 1/2 



x (s) — x (s — 1) 
a: (s + 1/2) - z (s - 1/2) 



(34) 



Thus, 



Vp (s) = a g " 1/2 [x (s) + x(s- 1)] + Pq~ 1/2 . 
Now, if one assumes that U is classic (see equation ([23])), or equivalently, 

V[p(s)W]=t(s)W, (35) 

one gets 

1 \ I , ... 1 



V\p(s)U] 



x (z) — x (s) 



-(s)U, 



x(z) — x (s) 



(36) 



According to Definition [31 one gets for the right hand-side of the equation (|36p 
the following relation 



t(s)U, 



1 



x (z) — x (s) 



U, 



is) 



x (z) — x (s) 



U. 



-js) - t(z) + t(z) 
x (z) — x (s) 



= r{z)S q (z) - bi (U, 1) - T(z)S q (z) - b lU q . (37) 
By Lemma [2] and taking into account the formula before ([3"2l , we have 



V S V, 



x (z) — x (s) 



V. 



q 1 ' 2 
ir(s,z) 
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Therefore, from Definition O for the left-hand side of equation (131)1 one gets 

1 \ i/9 I '. . P (s) 



V s [p(a)W], 



a; (z) — x (a) 



7r(s,z) 



Thus, 



p( s ) 

7r(s,z) 



p(s) -p(s - 1) + p(z - 1) 



n(s,z) 



= -q-^ 2 p(z-l)\7S q (z) + {U 



p(s)-p(z- 1) 

7r(s,z) 



Since, 



p (s) - p (z - 1) = a [x 2 (s) - x 2 (z - 1)] + /3 [a; (s) - x (z - 1)] 
From expression (|32[) we have 

p(s)— p(z — 1) a x(s) + x(z — 1) (3 1 



7r(s,z) 



Thus, 



p(s) -p(z - 1) 

7r(s,z) 



<7 x(s) — x(z) q x(s) — x(z) 



x(s) + x(z-l)\ ft 



q \ ' x (s) - x (z) 



-- (W, 



x (z) — x (s) 



a 



~[x{z)+x(z-l)] + ?-)S q (z). 
1 1/ 



Taking into account (|34[) the above expression transforms into the equation 
p(s) -p(z- 1)' 



-^= a(? - 1 U g- (? - 1 / 2 ^(z)Vp(z). 



7r(s,z) 

Therefore, the following expression for the left-hand side of ((36 



V [p (s) W] , g (z) * g (a) ) = P (* - 1) V5 9 (z) + S 9 (z) Vp (z) - aq- 1 ' 2 ^ 



= V\p{z)S q {z)]-aq- 1 ' 2 ul. 



(38) 



holds. 



Then, by using (|3S ]l . ([37]) . and (|38]l one gets 
V[p(z)S 9 (z)] = T{z)S q {z)+(aq- 1 ' 2 -b^ul 
= r(z)S q (z) + 



a 2 q 



-1/2 



-6 l5 - x (g - 1) - h 



and equation (|33|) holds. 

Notice that the statement holds since the converse implication, i.e. (|33|) = 
(|35[) . follows by the above chain of equalities but proceeding in reverse order. 
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Remark 2 Theorem]]] can be proved in a more elementary way, i.e., without 
using the assertions of Lemma [7] and Lemma fj| respectively. However, this 
procedure requires some cumbersome calculations. 



Below we highlight the main ideas of this procedure. Assume that IA verifies 
or cquivalently 



( V [(a(s) + t(s) V x( S | 1/2))U] , ) = (r(s)U, [s] q k) ) 



Hence, 



(U, [a(s) + r(s) V x(a + 1/2)] q 1 ^ 2 [k] q [s]^ + r(s) [,]<*> ) = 0. (39) 

From equation (|3"9")1 , after some straightforward computations, one gets a three- 
term recurrence relation involved the following g-moments u k _ 1 , u k , and u^.^. 
The resulting equation can be expressed in the following convenient way 



r fc + * fc = E k , 



(40) 



where 

r fe = - 



029 2 +-b 1 (q-l)q 2 



3fc-3 

-a 2 q 2 + 



aiq^ + ^hq^ 1 (q k 1 + q k - 1) + ^b q% l {q-l) 



a 2 qi- 1 [k - 1]J + (a x + \b iq k -'i) [k - l] q + (a + 



"fc-i' 



a 2 q k - 3 H q + l) + -b iq k - 2 (q 2 -l) 



4+i + { a 2 [<z* 2 (<? + i) [*],-g~ 



+ ai q * + ~h q 2 (1-q 2 ) [k] q + q~ 



1(9-1), , , 



and 



S fc = b iq k u q k+1 + (6igV [k] g + b )u q 



tedious calculations, one obtains 



Dividing T^, "J^, and Sj. by g fc [s]^ and summing from fc = to oo, after 



E 



fc>0 1 l S iq 



3 1 

«2g _i + -biq~ 2 {q - 1) 



+ [<r(a - 1) + r(s - 1) y - 1/2)] V5,(a), 



(41) 



fc>o 9 fc [*] 



a 2 <Z 5 (9 + l)+ 2 & i(!-9 2 ) 



+ V [a(s) + t(s) v x{s + 1/2)] S,(s), 



(42) 
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and _ 

5fc 



k>0 1 l b lq 



-hu q + T(s)S q (s). (43) 



Recall that from (|40|) one has 

r, ^ *, 



E _ fc r„i(fc+i) + E ' 



fe>0 y PJg fc>0 9 l S J g fe>0 1 l S \q 

Then, relations (|41l) - (|4"3")l yields the equation (f3"3"|) . The converse implication 
follows by the chain of obtained equalities but proceeding in reverse order. 

Observe that the more standard properties that characterize the g-classical 
orthogonal polynomials can be derived from a distributional difference equation 
(|35p . For this observation we refer to the approach given in [19] in the framework 
of a pure algebraic approach -used also in the above proof. Accordingly, from 
Theorem [T] follows that the g-classical orthogonal polynomials verify the hyper- 
geometric- type difference equation (|10[) . the three-term recurrence relation (|17[k 
their finite differences constitute an orthogonal polynomial family (|16|) and they 
can be expressed by a Rodrigues-type formula (|15[) . among other aforementioned 
properties. 



5 Examples 

In this Section, we will give an explicit expression -case by case- for the moment 
sequence ([^)) as well as for the solution of (|3"3"|) in terms of g-hypergeometric 
series, i.e. of expression (l24l) . 



5.1 g-Charlier case 

The q-Charlier polynomials are orthogonal with respect to the q-moment func- 
tional (I2TT) defined by the weight function (see e.g. [U eq. (87)]) 

u s 

e q [[l-q)fi] r g (s+l) 

where s 6 [0, oo), and e q (z) denotes the q -analogue of the exponential function 
(for details, see [XT]). 

Observe that the above function p(s) is a solution of equation (fT4"|) for the 
polynomial coefficients: t(s) = fiq 3 ^ 2 — q 1 ^ 2 x(s), a(s) = q s x(s), and 



cr(s) +t(s) xjx(s + i) = \iq 

These polynomial coefficients (see e.g. eq. (86)]) are involved in the charac- 
teristic equation (|33j) as well as the moment Uq. For determining this moment 
we now compute explicitly all the moments of (pH) for the above choice of p(s). 
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According to (p?!?|) one has 



g 1/2 v f s i (fe) = q 1/2 V M ' 

e q [(1 - ^ [ ^ T q (s + 1) e q [(1 - g) M ] ^ (q;g)s-k 



(1 - g) s - fc 



e 9 [(1 - ^ (g;g)n 
" (1-g)" 

Thus, the explicit expression for the fc-th moment associated with (f2Tj) is as 
follows 

_ M'e, [(1 ~ q)lA 
Uk ~ q ^e q [(l- q M> fc - U ' 1 '-" 

In particular 

q = e q [(1 - q)JI\ 
gV2 eg [(l-g) M ]- 

Finally, from (f24f one gets 



5 (z) = ji_ y (g; g) fe g Q [/*(! - g)g 1_ 



x ( z )fe> (g 1_z ;g) fe («;?)* 

since the g- falling factorial can be rewritten as follows (see (P2U10 

M <^> = (_!)*-£« (^; 9 ) fc9 ('- 1 )*-(S). (44) 
(1 - g) 

Hence, the g-Stieltjes function associated with the g-Charlier moment functional 
that solves equation (|33|) is given in terms of the following hypergeometric series 



S q( Z ) = — TT Ifl 

X (z) 



q 

,1-z 



q^(l-q)q 1 - z . (45) 



On the other hand, we know that (|4"5"|) must be equal to (f5U|) . This fact is easily 
established if one rewrites p(s) as 



and ([501) as 



(g;g) s ' q e q [(i 



^ (g 1 z ;g) s (g;g), 



C,«T 1/a ( g - ;9 ) s ( Q; ^ [(l-g) M9 ] s 



x ( z ) ^ (g 1 z ;g) s (g;g) s 
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Here we have used that F q (s + 1) = ppigji • 

Thus, the g-Stieltjes function associated with the g-Charlier moment func- 
tional also has the form 



S q (z) 



C q g-V * ( 
x(z) 



2fl 



q; (i - q)w 



Now, using the Heine's transformation formula (for details, see [12] page 16) 
a,0 



2fl 



q;z \ =e q [z\ npi 



q;az , a^O, \z\ < 1, 



under the assumption |(1 — g)/xg| < 1 one gets the desired equality: 



S q (z) = 



C q q-V * I 
x (z) 



■ 2<fl 



q; (i - q)m 



x(z) 



1<P1 



q 

_1— z 



- q)q x 



5.2 g-Kravchuk case 

The g-Kravchuk polynomials are orthogonal with respect to the q- moment func- 
tional (|2ip , where the orthogonalizing weight function 



p(s) = Q® r ( Jl —p s (l~ P ) N - s , 0< P <1, NeN, 

T q [s + l)T q (N - s + 1) 

is a solution of Pearson-type equation (|14p for the choice (see [1] p. 89) 

gl/2pg(g" - 1) g 1 / 2 (p(q - 1) + 1) 

r(sj = x(s) 



1-p 

cr(s) = (g- l)x 2 (s) + x(s). 



1-p 



By using the expression T q (s + 1) = g( 2 V 2 [s] 9 !, the above function p(s) can be 
rewritten as 

*>=^ ,,,M?,,2 wV' (i - p)N " 
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Hence, for the moments (|2"21 one gets 



M(")+i] y q s(N+D/2 mw W s(1 _ x 



AT- 



s — 



= 9 



TV 



s(AT+/c+l)/2_ 



s—k 



[s - k] q \[N - s] q !< 



^ h cr)-(;)-i][4VE 



^ g n(J V+fc+l)/2 [JV _ fc] g , n 



[n] g ![iV-fc-n]-l 
where we have used the following relations 

[iV] g ! = [7V]( fe) [iV - Jfc],!. 

Equivalently, 



P n (l-P) 



N—n—k 



N-k 



4= g i[»m-(M[jv](VE 



n(2fc+n+l)/2 



ri=0 



n 



p"(l-p) 



N-n-k 



where the g-binomial symbol is defined as (see e.g. [TTJ p. 24]) 



n 
j 



r,(» + i) 



r,(fc + i)r,(n-fc + i) 

Using the well-known relation (see O pp. 3-6]) 

n-l 



= q Hn-k)/2 



[k] q \[n-k} q V 



3=0 



n 
J 



J=0 



the above expression for the fc-th moment transforms into 

N-k-l 

^ = g i[2m-m-i] [iV ](v n (i-p+p? 3+fc+1 ) 



J'=0 



Finally 



where 



q _ "o 



(1 - q) fe ( _n_ 

p-1 



(1 / 



iV 
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Therefore, from (|24l) and f4"4")l one gets the following explicit expression for the 
g-Stieltjes function associated with the g-Kravchuk moment functional 

Sq[z) =a_y: {q ~ Nl q)k (<z; q)k { ~ 1)k 



pq 
p-i 



■ q 



(q;g)k 



This expression can be written in terms of hypergeometric series as follows 



S q (z) 



x (z) 



2</?2 



pq 



\ p-l 



pq 



N+l-z 



p-l 



(46) 



/ 



Previously, we have established a relationship between (|46j) and (|30|) . Now, 
aimed to check a similar relationship between (|46l) and (|30|) we rewrite p(s) as 



-JV. 



where 



P ( S ) = C q 

c q = (i- P ) 



(i-p)V^ /2 - 



(47) 



r 9 (iv + i) 

In (|47|) we have taken into account the following relation 
r g (iV- a + l) _ (-I)* (1 - q) S q&~ N 

r q (N + i) 

Thus, from one gets 



^ (z) = c ^[ 2 V (ff - ^ ^ (<r * ; *>• (^) f 



^ (q 1 z ;q)t 



{g;g)s 



Accordingly, the g-Stieltjes function associated with the g-Kravchuk moment 
functional also has the form 



S q {z) ^ q 



x (z) 



2^1 



I q- N ,q- z 


pq N+l 


V q 1 - 2 


J 



Now, using the Jackson's transformation formula (for details, see [H] p. 15) 
a, b 



2fl 



(az: o)„ 

(*;«)«, 



c, az 



q-bz , b^O, (48) 
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/ q- N ,q-* 


N+l ^ 






p- 1 

/ 



one obtains the aforementioned relationship 

C q q~^ (F 2 !^) 
N 



( q~ N ,q 



■ 2</?2 



/ q-\q 



■ 2V2 



x (z) 

5.3 g-Meixner case 



pq 



\ p—1 



pq 



N+l-z 



pq 
p-i 

\ 



pq 



N+l-z 



p-l 



P-I 



The g-Meixner polynomials are orthogonal with respect to the g-moment func- 
tional (|21|l defined by the weight function (see e.g. [H p. 314]) 



P(s) 



M*r g ( 7 + s) = y s (g 7 ;g), 
ry 7 )r ? ( s + i) Ji)m s] , (i- g )< 



•, < /i < 1, 7 > 0. (49) 



Recall that p(s) is a solution of Pearson-type difference equation for the poly- 
nomial coefficients (see e.g. [21 eq. (68)], up to the factor q — 1): 

t(s) = q 1 ' 2 ^ 1 - l)x(s) + t iq 2 ^[ 1 ] q , 
o{s) = (q- l)x(s) 2 + x(s). 

For the above orthogonalizing weight function and (|2"Tj) . let us compute the 
associated g-moments 



s>0 



n-him+A ^^(fc+2)/2-q) /2 m s (g 7 ;g) 



s>fe 



i[( fe + 1 )+i-*;(fc+2)] k K(fc+2)/2-("i fc V2 (g 7 ;g)»+fc 

^ ^ y [n],! (l-g)«+ fe ' 



n>0 



Considering the following elementary relations 

'n + k 



n(k + 2) 



n(5 - n)/2 



{q~'\q)n+k = {q~'\q)k{q 1+k \q)m 
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one gets 



fe-1/2 (g 7 ;g)fc „fc'Sp „n(5-n)/4 (g 7+fc jg)n 



fc-l/2 r ?(7 + fc) fcV^ n(5-n)/4 (g 7+fc j g)n 

r,(7) ^ „Vn W9I (i-g) n 



i 9 17 ' 1 n>0 

fnl„! = o~( 2 )^ 2 , the following relation 



holds. Thus, by the g-binomial theorem (see e.g. [TTJ eq. 1.3.2] ), 

u , = 1/2 (w 7+1 ;g) 00 r 9 ( 7 + fc) g 
fc (w;g)oc r,( 7 ) (/^ +1 ;g) fe ' 

In particular 

, _ x/2 (w 7+1 ;g) 00 

"o ^ g 



(Mg;g)oc 

In order to obtain an expression for the corresponding q-Stieltjes function asso- 
ciated with the q-Meixner moment functional one uses the following elementary 
relation 

r g ( 7 + fc) _ (q-y;q) k 



r,(7) (i-g) fe ' 

as well as (|2^|) and (El) , respectively. Thus, 

«§ ^ (g 7 ;g) fc (g;g) fc (-i) fc g© {m 1 -^ 

fe>0 



(50) 



x(z)^ n (g 1 — , 9 ) fe (Mg 7+1 ;g) fc (g;g) fc 



Equivalently, in terms of hypergeometric series 



x{z) 



g 7 ,g 



g;/V~ z ■ (51) 



Again, to establish an explicit relationship between (|5T|) and (|3U1) we rewrite 
(|49|) -based on equation (|50)) - as 

/ \ s (g 7 ;g) s 

p(s) = jU 



(g;g) s 

Thus, from (1501) one obtains the following g -Stieltjes function associated with 
the q-Meixner moment functional 

s r-) - 9 " 1/2 y (g 7 ;g) s (g _2 ;g)« (w) s 



^ (g 1 2 ;g)« (g;g) s ' 
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or equivalently, in terms of hypergeometric series 



S q (z) 



g-V 2 
x (z) 



2</ 3 l 



q 



Again, using the Jackson's transformation formula (1481) we have the following 
relation 



S q (z) 



9 



-1/2 



x(z) 



2<Pl 



q-w ( w 7+1 ;g) c 





( q\q 




2</?2 













q J ,q 



2</?2 



x (z) 

5.4 g-Hahn case 

The g-Hahn polynomials are orthogonal with respect to the g-moment functional 
(l2Tj) defined by the orthogonalizing weight in the interval [0, N — 1] (see [3j Table 
4.1]) 



(N + a-s) 



T q (s + l)T q (N-s) 



a,/3>-l, iVeN, 



where f(s) = q-^-^^-^^fis; q) if < q < 1, or f(s) = /(s;^ 1 ) if q > 1, 
being f(s;q) = (1 - q) 1 -* I\ k > (l - q k+1 )/ U k >o( l ~ 9 s+k )- The above function 
p(s) is a solution of equation (fT4|) for the polynomial coefficients (see [3]): 

a(s) = - (Z -( Ar+Q )/ 2 x 2 ( s ) + q-^ 2 [N + a] ? x(s), 

r(s) = -q-^ +2 - N ^ 2 [a + (5 + 2] fl i(«) + g"+^ +1 [/3 + 1],[JV - 1],. 

For calculating the q-moments f|22[) we use 



M fe = q 



= q 



v „(h±|±2) s r ,(*) T q ( s + /3 + l)T q (N + 
^ q [lq \8] q l[N-l-8l 



s—k 



a — s) 
J 



i] ^ (2±£±fc±2) s r g (s + /3 + l)T q (N + a — s) 
^/ 2 [ S -fe] ? ![JV-l- s ],! 



Hence, 



JV-fc-l 



,.8 _ n x(k,a,/3) 

a k — q 



" n( - Q+/3 + fc+2 - ) T g (/3 + fc + 1 + »)r g (A 
^ 9 2 [nl,![JV-*-l-»j,, 



n=0 



(JV + a - fc - n) 
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where X (k,a,f3) = k (h±£±*±E) _ i [ (fe +i } 

.ft _ g 



, or equivalently 



[N-k-iy 

N-k-1 

q iP(a,/3,n,N) 

n=Q 



N-k-1 
n 



T q (fi + k + 1 + n)T q (N + a-k-n). 



For brevity we have introduced the notation: ip (a, j3, n, N, k) — n(a + ft + 2k + 
3 + n-N)/2. 

Taking into account 



(g/»+*+l;g) 



(1-?)' 



Tqjfi + k + 1 + n) _ ( ( 3 + fc +n )( ( 3 + fc +n _i)/ 4 

r, {fi + k + 1) 

r g (A^ — fc — n + a) _ (JV-fc-l-n+a)(W-fc-n+a-2)/4 

r g (a + 1) 
the g-moments becomes 

„ _ r 9 (a + 1) r q (13 + fc + 1) g^w) 



(1 _ ^JV-fc-l-n ' 



jiV— fc— I 
iV-fc-1 

E * n(a 

n=0 



+ 1) 



N-k-1 
n 



^^N-k-l-A^^) 



where 

2(t>(a, /3, AT, fc) + 1) = a (2k - N + 1) + N (fc + 1) - 
Using the following well-known relation (see [TTJ p. 25]) 



fc\ /JV\ /or\ A3 



E 



n 



one gets 



,9 - 



v(a,P,N,k) 



T q (a + 1) (/? + fc + 1) g 



' N-k-1 



In particular 



_ r g (a + i)r g G9 + i) g "(°^°) , +/3+2 . , 
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On the other hand, using relations 
lV/3 + fc + l) = (q fj+1 ;q) k 

r,C8 + i) (l-q) k 

r^v-fc) = _ w _ fc _ lKiV _ fc _ 2)/4 (-i) fc (i - g ) fc gft-wffl 

T q (TV) = g (^-D(w-2)/4 
we obtain 

M 9 - u « *< w+a HS) (^ +1;g )fc( gl ' W;g )fc 

* °* ( } (i-s) fe (^+W 

Finally, using ([M]) and (HH) one gets the g -Sticltjcs function associated with 
the g-Hahn moment functional 



k>0 



Q ( z ) = "0 V 

?M X(Z)^ (gl-« ;g)fc (g«+fl+2 ;(?)fc (?;(?) 



or equivalently -in terms of hyper geometric series- 



u 

Sq(z) = —T^3<P2 



\ q l-z a+/3+2 



q; q N + a - z 



Now, aimed to find an equivalent relation for the above g-Stieltjes function one 
uses the relations 

r 9 (s + i) = ^^h, 
I 1 - q) 

and 

r q (s + p + i) (s+ms+0 _ iyA (/ +1 ;g) s 
M/3 + i) 9 (i- q ) s ' 

as well as 

T q (N + a-s) = - (N+a -s-i)(N +a -s-2),4 (~l) a (1 - g) S g (1 - JV - a)a+ (') 

r g (7V + «) 9 (q 1 - N - a ;q) s 

Thus, we can rewrite p(s) as follows 

PU q T q {N) tf-^q^^qh' 

where 2v(a, j3, N) = a (1 - N) - (°) - Then, from flU) one gets 

5 ,. V - v (^ +1 ; g ).(g 1 - y ; g ) a ( g -'; g ). g' Ar + a . N 
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where 



Therefore 



Cq - Q T q (N) 



S q {z) = 

x (z) 



3^2 



g l-JV-a 



g;g 



(52) 



Finally, taking into account the relation 



°" (g^g)^ 

and using the transformation formula (see Theorem 12.4.2 

. 6" (d/6;g)„ 



3^2 



q n ,a,b 
c, <i 



g;g 



■ 3f2 



q n ,b,c/a 
c,g 1_n 6/d 



taking n = N— l,a=q z , b = q l3+1 , c = q 1 z , and c? = q 
transforms into the following expression 



l-N—a 



relation J52 



S q (z) 



x (z) 



3<P2 



q? + \q X - N ,q 



n l-z n a+i3+2 



g;g 



N+a- 



where relation (aq- n ; g)„ = (g/a; g) n (-a) n g -™(«+i)/2 ( see [7J p. 304, (12.2.10)]) 
have been used. 



6 Conclusions and future directions 

Based on the accumulation of analytic and algebraic properties that characterize 
the classical discrete orthogonal polynomials -often seemed to be unrelated- the 
need for structure and classification of them has constantly arisen as a central 
question in the orthogonal polynomial theory [9j [15l [TH [20] . In this paper we 
address this question proving, in Theorem [lj that a sequence of ^-polynomials 
orthogonal with respect to g-moment functional (|2"Tj) is classical iff the associated 
Stieltjes function given in terms of q -falling factorial basis satisfies a first order 
non-homogeneous q -difference equation; the proof is given in a constructive 
way using a theoretical background based on the theory of linear functional 
deeply studied in [16] and [19]. We show, in Theorem Q] that the verification 
of the aforementioned difference equation constitutes a new characterization 
of a discrete orthogonal polynomials on the non- uniform lattice x(s) = (q s — 
1) / (q — 1). However, more general situations demand special attention. For 
instance, the g-orthogonal polynomials on the lattice x(s) = c\q s + ciq~ s + C3, 
(q € R + \{1}), where the constants Ci £ K (i = 1, 2, 3) are constants independent 
of s, must be considered in the same fashion as here. This paper outlines the 
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important points and techniques to be followed in such investigations aimed to 
characterize those families of g-polynomials. 

Finally, more general systems of g-moment functionals must be analyzed. In 
this direction, the g-semiclassical orthogonal polynomials could be an interesting 
challenge to be considered; in particular when Dirac masses are added to q- 
moment functional. An analogous result to those given in Theorem Q] played a 
crucial role in the computation of the class of the semiclassical functionals given 
by a perturbation via the addition of Dirac masses (see [TU]). 

In closing, to the best of our knowledge, there is not in the literature any 
explicit expression for the associated Stieltjes functions in the q -falling factorial 
basis given in terms of hypergeometric functions. 
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